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Abstract 

We study the longitudinal polarization of the tau lepton in B — )• Dri) decay. After discussing 

<^ possible sensitivities of r decay modes to the r polarization, we examine the effect of charged Higgs 

"O boson on the r polarization in i? — )• Dtv'. We find a relation between the decay rate and the r 

(N 

i— H polarization, and clarify the role of the r polarization measurement in the search for the charged 
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I. INTRODUCTION 

The quarks and leptons in the third generation are important clues to new physics beyond 
the standard model (SM). Because of their larger masses their couplings to the electroweak 
symmetry breaking (EWSB) sector are relatively strong, and thus their interactions are 
potentially sensitive to new physics that modifies the EWSB sector of the SM. 

The minimal supersymmetric standard model (MSSM) is well-motivated and attractive 
among several candidates of such new physics. Its Higgs sector at the tree level corresponds 
to the two-Higgs-doublet model (2HDM) of type II and contains a pair of charged Higgs 
bosons H^ in the physical spectrum. In the 2HDM of type II, the first Higgs doublet 
couples to the down-type quarks and the charged leptons, while the second one does to the 
up-type quarks. Then, the interaction of the charged Higgs boson with fermions contains 
terms proportional to r/i/tan/?, where irif denotes a down-type quark mass or a charged 
lepton one, and tan/3 = ^2/^1 with fi(2) being the vacuum expectation value of the first 
(second) Higgs doublet. Consequently, the charged Higgs contributions to the amplitudes 
of tauonic B decays involve terms proportional to m;,m^ tan^/3, and is enhanced if tan/3 
is large. Several theoretical and experimental studies on tauonic B decays have been done 
motivated by this observation. 

The branching fraction of the pure tauonic B decay, B — t- rz/, is measured as (1.7±0.6) x 
10"^ (BABAR) [1] and (1.65l[]J^+[]:^^) x 10"^ (Belle) ^2]. Combing them, the Heavy Flavor 
Averaging Group (HFAG) obtains (1.67 ±0.39) x 10"^ (HFAG) [3]. 

The theoretical estimation of the decay rate of -B — > rP including the charged Higgs 
effect is straightforward |1], but suffers from significant uncertainties in the uh element of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix OE] \Vub\ and the B meson decay constant fs- 
Taking the ratio oi B ^ ri? to B ^ fiu does not help, since the lepton universality in the pure 
leptonic B decays is not spoiled by the charged Higgs effect [4]. Using the branching fraction 
given by the HFAG, iKbl = (3.95 ± 0.35) x lO'^ [7] and fn = (190 ± 13)MeV |B], we obtain 
the allowed region of 95% C.L. for the charged Higgs parameter as tan (3 /rriHi < 0.11 GeV~^ 
and 0.24 GeV^^ < tan/3/mj:/± < 0.31 GeV~^, where m^i denotes the mass of charged Higgs 
boson. 

Semitauonic B decays are more complicated than the pure tauonic B decay. However, 
there are several observables in them besides branching fractions, e.g. decay distributions 



and r polarizations. This is one of the reasons that we study B — )■ Dtu in the present work. 
Another reason is that it is known to be the most sensitive to the charged Higgs among 
several semitauonic B decays studied so far [H]. 

The charged Higgs effects on the branching fraction, the q^ distribution, and r polar- 
izations in i? — 7- Dri? are investigated theoretically in the literature [9H2T]. The present 
experimental results on the branching fraction of i? — ;■ Dtu are given by BABAR and Belle 
collaborations: 

^ ~^ ^ '^^'' -0.416 ±0.117 ±0.052 BABAR [22] , (1) 



B{B^ Di-Ui) 



and 



B{B -^ D T V) ^ .0+0.22+0.06 R„ll„l9ql (n\ 

B{B^ ^ D-i^u) - 0-48-0.19-0.05 Belle [23] , (2) 

B{B+ -^ D^T+v) = (0.77 ± 0.22 ± 0.12)% Belle [21] . (3) 

Averaging them, we obtain the branching-fraction ratio R as 

where B{B^ -^ D^C^v) = (2.15 ± 0.22)% is used [23]. A more precise measurement with a 
few percent error in the branching fraction is expected in a super B factory |25j . 

The daughter r's are identified by successive r decays in the experiments: r — )■ iuu 
[i = e, /i) is used in Ref. [22| and Ref. [231 while both r — > iuu and r — > ttz/ are used in 
Ref. [211 The distribution of r decay products in i? — )■ Dri? is also sensitive to the charged 
Higgs boson [20j . It is illustrated in Ref. [201 that the distribution of the angle between the 
momenta of the D meson and the pion in r — t- ttz/ in the B rest frame depends on the 
magnitude and the complex phase of the charged Higgs coupling. It is also expected that 
the information on r polarizations, which are affected by the charged Higgs boson [91 IT111T3] . 
is encoded in the decay distribution of successive r decays. 

Effects of r polarization on the r decay distribution are well-studied for many r produc- 
tion processes, such as e~e+ — > t^t^ [2SHS0] . Z° — )■ t^t^ |SIH31], Higgs decays [Ml ES], 
other heavy particle decays [36VI38] . and the i/T--nucleon scattering [39]. It is shown in the 
literature that we can decode r polarizations from appropriate r decay distributions in these 
processes. We apply a similar method to B — t- Dtu in the present work. 



In this paper, we study the r longitudinal polarization in B — ^ Dtv and clarify its role 
in new physics search with the main interest in the charged Higgs boson. It turns out that 
the r longitudinal polarization combined with the branching fraction gives us a valuable 
hint for new physics. In Sec. [IT], we examine possible sensitivities of r — ;■ ttz/ and r — )■ Ivu 
to the r polarization in i? — )■ Dtv. Then, we summarize the charged Higgs effects on the 



decay rate and the r polarization in Sec. Ill We show our numerical results including a 



relation between the decay rate and the r polarization in Sec. IV This relation results from 



a distinctive nature of the charged Higgs interaction. Sec. V is devoted to conclusions. 



II. TAU POLARIZATION AND ITS DECAY DISTRIBUTION 

In this section, we illustrate how to extract the r polarization in B ^ Dtv using suc- 
cessive r decays and examine possible sensitivities in experiments. It is possible to define 
two distinct and independent r polarizations, namely the transverse polarization and the 
longitudinal one. The transverse polarization is known to be generated by T violating inter- 
actions and/or final state interactions [TitiTS] . Both interactions are small in the SM. While 
T violating effects induced at one- loop level may be sizable in the MSSM. The longitudinal 
polarization is supposed to be sensitive to the chiral structure of the relevant interactions. 
The interaction of the charged Higgs boson, which is our main concern in the present work, 
has a different chiral structure from that of the W boson. We concentrate on the longitudinal 
polarization in the following. 

The r longitudinal polarization depends on the frame in which it is defined. We employ 
the frame in which the spacial components of the momentum transfer q^ = p'^ — p'^ vanish, 
where p'^ and p^ are the four-momenta of the parent B meson and the daughter D meson 
respectively. We refer to this frame as the q rest frame. Note that the q rest frame is 
accessible in the e^e~ B factories provided that the tag-side B meson is fully reconstructed. 
Incidentally, the q rest frame corresponds to the center of mass system in e~e+ — > t^t'^ in 
the sense that both are the center of mass system of the lepton and anti-lepton pair. In this 
way, the choice of the q rest frame turns out to be reasonable. 

We use a coordinate system in the q rest frame such that the direction of the B and D 
momenta is the z axis, and the r momentum lies in the x-z plane. Then, we parameterize 
the r momentum as p^ = Er{l, PrSinOr^O, l3rCos6r), where Er = {q^ + m1)/{2^/<f), j3r = 



a/1 — rri^/E^ . The helicity amplitude of 5 — )> Dru is denoted as A^'^^(g^, cos^t-), where 
A^ = ± designates the r hehcity defined in the q rest frame and the neutrino hehcity 
is assumed to be negative. The exphcit form of A^'^^ is given in the next section. The 
differential decay rate of -B — ?> Dri? for a given r helicity At is written as 

dTx^ = J—\M''-{q\coser)\^d<!>s, (5) 

ZniB 

where the three-body phase space d^^ is given by 

rf$3 = ^^^ f 1 - "4) dq'd COS 9^ , (6) 

and Q± = {rriB ± itld)'^ — q^. The r longitudinal polarization in the q rest frame is defined 

as 

dr\ fdV, dV, 



where dV /dq^ = dVj^/dq^ + dT^/dq"^, and we integrate over cos^^ since it is difficult to 
determine the direction of the r momentum at the B factories in contrast to the case of 
Z^ — 7> r^r^ at LEP and SLC. Furthermore, we introduce the average r polarization, 

where F = r_|_ + r_ is the decay rate of -B — ?> Dru. Though the average polarization holds 
less information, it is still sensitive to the charged Higgs as we will see below and supposed 
to be useful for experiments with limited statistics. 

Measuring the r polarization in addition to the decay rate summed over the r helicity, 
we can separately extract the decay rates of ffxed r helicity. This gives us a clue to new 
physics as we will see later. 

The r in i? — 7> Dri? is identiffed by r — )■ ttz/ or r — > iuu {i = e, /i) in the present B factory 
experiments as mentioned above. Accordingly, we see how these r decay modes work as r 
polarization analyzers. The differential decay rate of the decay chain B — )> Dri? followed by 
r — 7> TTZ/ (r — )■ iuu) is written as 

|f-B.,4l/tf.C) + P.tf)«tf.C)l (9) 

where ( = ET,(^e)/Er with Et,(^i) being the 7r(£) energy in the q rest frame and 3^,(^1) denotes 
the branching fraction of r — )■ ttz/ (r — )■ 



The functions / and g for r — t- ttz/ are well-known and given by 



/(g^C) = l//3., ^(g^C) = (2C-l)//3^ 



(10) 



where we neglect the pion mass for simplicity, and the range of C is (1 ~ /3t)/2 < C ^ 
(l + /3r)/2. 

As for r — )■ ^z/i/, ignoring the E mass, the decay distribution is described by 

16 C' 



/(?^c) 

9{q\0 



3 (1 - /3^)3 
16 C' 



3 (1-/3^)3 



[9(1 - Pi) - 4(3 + Pl)(\ 
/3.[3(l-/32)-16C], 



(11) 
(12) 



for < C < (1 - /3r)/2, and 

1 + /3r - 2C 



/(?',C) 
^(g',C) 



3/3,(1 + /3.)3 
1 + /3, -2C 1 



3/3,(1 + /3,)3/3, 



[5(1 + /3,)2 + 10(1 + /3,)C-16C'], 

[(1 + Prf + 2(1 + /3,)C - 8(1 + 3/3,)C'] 



(13) 
(14) 



for (1 — /3r)/2 < C ^ (1 + I^t)/'^- Eqs. (13) and (14) reduce to the more familiar functions 
in the coUinear limit /3t- — > 1, see e.g. Ref. |3H 

We can determine PiiQ'^) by measuring the C distribution for fixed g^ in Eq. (pi). The 
statistical uncertainty of the ideal experiment is given by 

1 



SPiiq' 



(15) 



where N(q'^) is the number of signal events for fixed q^ (or in a bin of g^, more practically) 

1/2 



and 



S{q' 



dC- 



9\q\0 



f{q'X) + PL{q')9{q',0 
For the average polarization Pl in Eq. (p|, we obtain 



(16) 



SPr 



NS 



(17) 



where N is the total number of signal events and the average sensitivity S is given by 

.-1/2 



S 



dT 



dq'^,S-\q^) 



(18) 



Assuming the SM and neglecting the uncertainties in the form factors discussed in Sec. IV 



we obtain S = 0.60 and 0.23 for r — )■ ttz/ and r — )■ iuu respectively. These values vary 



less than 20% even in the presence of charged Higgs boson taking the constraint from the 
branching fraction into account. 

The expected uncertainty in P^ is 5Pl ~ 0.4 with A^ ~ 100 for r — )■ Cvv^ which corre- 



sponds to the present experimental status [22], |23]. As for r — )■ ttz/, 6Pl ~ O.S/A/eTr/e^^, 
is expected in the present experiments, where the branching fractions of r — )■ ttz/ and 
r — ;■ ivv are taken into account, and e.,^(t) represents the efficiency of the r — ;■ ttz/ (r — > Cvv) 
mode. At the super B factory with integrated luminosity of 50ab~^, A^ ~ 2000(3000) for 
r — 7- TTZ/ (r — )■ Cvv) is obtained based on the Monte Carlo simulation in Ref. 123 and thus 
5Pl ~ 0.04(0.08) is expected^ 

III. HELICITY AMPLITUDES AND DECAY RATES 

In the presence of charged Higgs boson, both the W boson and the charged Higgs boson 
contribute to the helicity amplitude of 5 — t- Dtv. We describe their contributions in turn. 
The W boson exchange amplitude A^^ is written as |10l HI] 

MW.^osOr) = ^K,5^r/,^i/A,,^t ' (19) 

Xw 

where Vcb is the cb element of the CKM matrix, \w = ±, 0, s denotes the virtual W helicity, 
and the metric factor rjx^ is given by ri±^o = 1 and rjs = —1- The hadronic amplitude Hx^ 
that represents the process B — ^ DW* is defined by 

HxM) = e*^i\w){DipD)\cril - l,)h\B{pB)) , (20) 

where e^{Xw) is the polarization vector of the virtual W boson. The leptonic amplitude 
L^^, that represents the process W* — )■ tv^- is defined by 

LllXq\coser) = e,{\w){r{prAr)i^r{p.)\fY{l-l,)vM. (21) 

Here, we introduce the hadronic form factors h±{w) |42j . 

{D{v')\cYb\B{v)) = y/niBniD [h+{w){v + v')" + h_{w){v - vY] , (22) 

^ We assume that efficiencies of r — > ttz^ and t ^ pv are the same. 



where v^ = p'^/niB, v'^ = p'^/niB and w = v ■ v'. The hadronic amphtudes are written in 
terms of these form factors: 

iJ±(g2) = 0, (23) 

Ho{q') = .M^^ —-^^'' ^ V^^^^^IViH, (24) 

V 1 — ^rw + r^ 

HM^) = ,/^^^^^ -=^L=={w + 1) Si{w) , (25) 

V 1 — 2rw + r^ 

where r = m£)/mB-, and 

V,{w) = h4w)-^h4w), (26) 

S,{w) = h^{w)-\±^^'^h.{w). (27) 

In the heavy quark hmit (HQL), h-^{w) reduces to the universal form factor known as the 
Isgur-Wise function C,{w) with the normahzation ^(1) = 1, and h_{w) vanishes [33j. Vi{w) 
and Si{w) also reduce to the Isgur-Wise function in the HQL. 
The required leptonic amplitudes are explicitly given as 



Lo(g2,cos^,) = -2y^q^y^l-ml/q^smer, (28) 

L+{q'^, cos Or) = 2m^v^l-m2/g2 cos^^, (29) 

L7(g^cos0,) = 0, (30) 

Lfiq^cose,) = -2m,Vl-m2/g2. (31) 



Note that the leptonic amplitudes other than Lq disappear for massless leptons, and thus 
the form factor that appears in i? — )■ Dii? {£ = e, fi) is only Vi. 

The helicity amplitude of the charged Higgs exchange is written as [9] 

M'^^ (g^ cos e^) = ^K, ^t J Hn L'^ , (32) 

where the hadronic amplitude Hr is defined as 

Hniq') = {D{pD)\ca + l^Wips)) , (33) 

and the leptonic amphtude is 

L^^iq^COsOr) = {T{pr,Xr)MPu)\f{^~l5K\0) . (34) 



The model-dependent coupling factor tp is given as t^ 
while 



ti 



tan^ 13 



tan^ /3 in the 2HDM of type II, 



(35) 



^ (l + £:otan/3)(l + £:^tan/3) 
in the MSSM, where Eq and Er represent radiative corrections pTl H^ . Using the equations 
of motion, we relate the hadronic and the leptonic amplitudes of the charged Higgs exchange 
to those of the W exchange with A^y = s as 



nibHR 



H. 



rrirL 



K 



72 L^ 



(36) 



1 -r„. 

where r^ = rric/mi,. Note that the charged Higgs contributes only to the amplitude of 
A^ = + and changes the r longitudinal polarization as well as the branching fraction. 

Substituting the total helicity amplitude Ai^^ = Ai^^+Ai]^ into Eq. ([5| and integrating 
over cos 6'^, we obtain 

dVx^ Gp\Vch\ 



dq^ 



128TT^m% 



Vq^- 






FxAq' 



where 



F-iQ') 
F+{<f) 



\q'\Mci')\\ 



,2\|2 



m^ 



\H,{q')V + 



^p q' 



m 



H 



2m2 



\Hs{q')\ 



(37) 

(38) 
(39) 



We clearly see the negative interference between the charged Higgs contribution and the 
standard W boson one as far as ti is positive. Once the form factors Vi and Si are given. 



we can evaluate the decay rate T\^ by integrating Eq. (37) over q . 



As noted above, the charged Higgs boson contributes to the rate of At = +, not to 
\r = —■ We can test this peculiar feature of the charged Higgs boson by measuring both 
the spin-summed decay rate and the r longitudinal polarization. 

IV. NUMERICAL RESULTS 
A. Form factors 



We employ the following ansatz for Vi{w) 



Vi{w) = Vi(l) [1 - Spjz + {5l.pl - 10.)z^ - {252.pl - 8A.)z^] 



(40) 



where z = (y/w + 1 — V2)/{y/w + 1 + a/2). Since Vi{w) governs B — t- Diu as we mentioned 
above, the slope parameter pf is determined by the experimental data of the q^ distribution 
in B ^ Diu. The recent analysis by the HFAG gives pj = 1.18 ± 0.04 ± 0.04 [3j. To 
be conservative about the uncertainties in the form factors, we combine the above errors 
linearly in our numerical work. 
We parameterize Si{w) as 

SiH = [i + AH]viH, (41) 

where A{w) denotes the QCD and l/wig corrections. We estimate the next leading order 
QCD correction following Ref. 161. In the numerical calculation of the QCD correction, we 
use rrib = 4.91 GeV and rric = 1.77 GeV as the pole masses of the bottom and charm quarks 
respectively |i47j, and as{mz) = 0.118 for the running strong coupling [7]. As for the l/rnq 
corrections, we take them from Ref.|45jand use A = 0.48 GeV for the mass difference between 
a heavy meson and its constituent heavy quark. 

Since the MS scheme is employed in the calculation of the QCD corrections, we should 



use the same scheme for the quark masses in Eq. (36). Thus, using the MS masses m;,^c(/^), 
Vm = rhc{p)/rhf,{p). Note that r^ is independent of the renormalization scale p as it should 
be. We use r^ = 0.21 in the following numerical calculations [47J. 

The analytic formula of A{w) is rather cumbersome and a detailed discussion on it is 
beyond the scope of this work. We only present an approximate expression 

AH = -0.019 + 0.041(w - 1) - 0.015(u7 - 1)^ , (42) 

which is as good as 3% in the physical range of w. In the following numerical results, we 
assume ±100% error in the estimation of A{w), that is, we replace A{w) by aA{w) and vary 
the uncertainty factor a from to 2. 

B. Decay rate 

Although the effect of charged Higgs on the decay rate is well studied in the literature, we 
present our numerical result to summarize the present status. It is convenient to introduce 
a normalized decay rate for each value of A^, 

^A. = ^ , (43) 

10 



ct^ 




0.2 0.3 0.4 

t^lmu± (GeV-^) 



FIG. 1: The branching-fraction ratio R = B{B — )■ DT^Ur)/B{B — )■ DH^ui) as a function oitp/mu± . 
The dark shaded (magenta) band represents the theoretical prediction including the uncertainties 
due to Pi and a. The light shaded (light blue) horizontal regions show the present experimental 
bounds at la and 2a. 

where F^ = F_|m^=o is the decay rate of -B — )■ Diu. We expect that several uncertainties 
(both theoretical and experimental) tend to cancel by taking the ratio of the decay rates. 
In particular, the uncertainty in |K;fe|Vi(l) disappears in the theoretical calculation. The 
branching-fraction ratio defined in Eq. ^ is given by i? = -R+ + R-. 

In Fig. fll we show the branching-fraction ratio i? as a function of ti^/niui, the control 
parameter of the charged Higgs effect. Hereafter, we take t^ to be real and positive. The 
dark shaded (magenta) band represents the theoretical prediction with the uncertainties in 
pf and a. The present experimental bounds corresponding to Eq. (E| are also shown in the 
figure by the light shaded (light blue) horizontal regions. A few comments are in order: 

1. The SM prediction is -R|sm = 0.302 ±0.015, which does not contradict with those in 
the literature [T8tl20]. 

2. The present experimental result is consistent with the SM, but it seems slightly larger 
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FIG. 2: The r longitudinal polarization Pl as a function oftf3/mjj±. The narrow shaded (magenta) 
band represents the theoretical prediction with the uncertainties due to pf and a. 

than the SM prediction. 

3. The allowed regions of tp/'mH± are given as tp/mH± < 0.14 GeV^^ and 0.42 GeV^^ < 
tp/^nnH± < 0.50 GeV~^. The latter region, in which the charged Higgs contribution 
dominates over the W boson contribution, is practically excluded if combined with 

B^ — )■ TV. 

C. Polarization 

In Fig. [2| the r longitudinal polarization in the q rest frame is presented as a function 
of tp/mH±. The width of the band shows the uncertainty in the theoretical calculation 
corresponding to p\ and a. The SM prediction turns out to be Pl = 0.325 ± 0.009. The 
theoretical uncertainty is remarkably small and dominated by the a factor. The expected 
statistical uncertainty in the super B factory is 5Pl ~ 0.04 and larger than the uncertainty 
in the SM prediction. 
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FIG. 3: The relation between the branching- fraction ratio R and the r longitudinal polarization 



Pl. The dark shaded (magenta) band represents the relation in Eq. (44) with the error. The 
present experimental bounds on R at la and 2cr are also shown by the light shaded (light blue) 
horizontal regions. The (blue) crosses indicate the theoretical predictions on R and Pl for several 
values of tji/mij± including the SM. 

D. Relation between R and Pi 

The decay rate and the r longitudinal polarization are independent observables in general. 



However, as mentioned in the last paragraph of Sec. Ill, they are related in the case of the 
charged Higgs because of the specific chiral structure of its interaction. It is straightforward 
to find 

/?(1-P^) = 2i?_ = 0.204 ±0.008, (44) 



where i?_ is determined only by the standard W boson contribution as seen in Eq. (38) 



We present this relation in Fig. p^ as the dark shaded (magenta) band with the error. The 
light shaded (light blue) horizontal regions show the present experimental bounds on R at 
la and 2a. The present experimental result on R implies 0.15 < Pl < 0.64. The theoretical 
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predictions on R and Pl for several values of tfs/mH± including the SM {tp/mH± = 0) are 
also indicated by the (blue) crosses. The leftmost cross is the turning point regarding the 
curve as a trajectory parameterized by tp/mH±. Incidentally, the two-fold ambiguity in 
tp/fT^H^ apparently remains. But, it can be solved combining with B — )■ tv. 



Eq. (44) provides a crucial test for the charged Higgs ansatz. If a set of R and Pl is found 
out of the dark shaded (magenta) band in Fig. Isl it immediately signifies the existence of 
new physics other than the charged Higgs. On the other hand, if one finds it within the 
band, but away from the SM prediction, it means that the new physics contributes to r_|_ 
and not to r_, and strongly suggests the charged Higgs. 

V. CONCLUSIONS 

We have studied the r longitudinal polarization in the q rest frame in i? — t- Dtv. The 
r polarization is measured through the distribution of subsequent r decays. The expected 
statistical uncertainty at the super B factory is 5Pl ~ 0.04(0.08) for r — t- ttz/ (r — ?■ Ivv). 

Then, we have examined the effects of the charged Higgs boson to the decay rate and 
the r polarization in B ^ Dtv. It turns out that the allowed ranges of the charged Higgs 
parameter for the present value of the branching fraction are tp/mH± < 0.14 GeV~^ and 
0.42 GeV~^ < tjs/mjj± < 0.50 GeV^^, and the uncertainty in the theoretical calculation of 
the r polarizations is notably small. 

Furthermore, we have found that the r longitudinal polarization Pl is uniquely related 
to the branching-fraction ratio R in the presence of the charged Higgs effects. This relation 
reflects the specific feature of the charged Higgs interaction. The present experimental result 
R = 0.40 ± 0.08 implies 0.15 < Pl < 0.64. If a deviation from the SM is found in R, the r 
longitudinal polarization will provide us an important information on the new physics. 
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